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Annihilation ol dark matter particles in cosmological halos (including a halo of the Milky 
Way) contributes to the diffuse gamma-ray background (DGRB). As this contribution will appear 
anisotropic in the sky, one can use the angular power spectrum of anisotropics in DGRB to constrain 
properties of dark matter particles. By comparing the updated analytic model of the angular power 
spectrum of DGRB from dark matter annihilation with the power spectrum recently measured from 
the 22-month data of Fermi Large Area Telescope (LAT), we place upper limits on the annihilation 
cross section of dark matter particles as a function of dark matter masses. We find that the current 
data exclude {av) > 10"'^^ cm^ s~^ for annihilation into 66 at the dark matter mass of 10 GeV, 
which is a factor of three times larger than the canonical cross section. The limits are weaker for 
larger dark matter masses. The limits can be improved further with more Fermi-LAT data as well 
as by using the power spectrum at lower multipoles [l < 150), which are currently not used due to 
a potential Galactic foreground contamination. 

PACS numbers: 95.35. +d, 95.85.Pw, 98.70.Vc 



I. INTRODUCTION 

Understanding the identity and nature of dark mat- 
ter, which makes up more than 80% of the total mat- 
ter density in the Universe, is a major goal of modern 
physics and cosmology. The most promising candidate 
for dark matter is the weakly interacting massive par- 
ticles (WIMPs), with which one can naturally explain 
the observed dark matter density using a simple thermal 
freeze-out argument T . If dark matter particles annihi- 
late into standard model particles, as expected for most 
WIMP scenarios, one can indirectly detect and constrain 
properties of dark matter particles [I]-^. In this paper, 
we shall focus on high-energy (1-50 GeV) gamma-ray 
photons produced by the cascade of annihilation prod- 
ucts. 

Dark matter annihilation occurs in all cosmologi- 
cal halos including a halo of the Milky Way, and 
thus contributes to the diffuse gamma-ray background 
(DGRB) .4r8 ■ Due to the large-scale structure of 
the Universe, the observed gamma-ray emission appears 
anisotropic in the sky in a predictable manner, which 
makes it easy to identify the dark matter origin of high- 
energy gamma rays in the sky [5] (also see Refs. [TDHIS])- 

Recently, the Fermi-LAT collaboration has measured 
the power spectrum of DGRB anisotropy from the 22- 
month of data [27j. They have detected significant ex- 
cesses of the angular power spectrum over the shot noise 



of photons for a multipole range between £ = 155 and 504 
and for multiple energy bins. A further study shows that 
most of these excesses come from unresolved blazars [35] . 
Subtracting the estimate of the blazar contribution, we 
have upper bounds on the residual anisotropy of DGRB. 

In this paper, we use the upper bounds on the power 
spectrum to constrain the annihilation cross section of 
dark matter particles. For this purpose, we update our 
theoretical framework for computing the angular power 
spectrum presented in Refs. [SJ [TUl HO] as follows: 

1. We use the results from recent numerical simula- 
tions (e.g., Ref. [IS]) to model the mass function 
and spatial distribution of subhalos within a given 
host halo. 

2. We include contributions from both the extragalac- 
tic dark matter halos and the Galactic dark matter 
subhalos. 

3. We also include the cross correlation between dark 
matter annihilation signals and blazars. Although 
this term was often ignored in the literature (except 
for Ref. [IQ]), one should include this term for self- 
consistency: the same halo hosting a blazar also 
contains annihilating dark matter particles, and 
there is a spatial correlation between halos hosting 
blazars and those not hosting blazars but contain- 
ing annihilating dark matter particles. 
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This paper is organized as follows. In Sec. [TTj we 
present the predicted angular power spectrum of DGRB 
from extragalactic dark matter halos, and compare this 
to the data to find constraints on the annihilation cross 
section. In Sec. |III[ we discuss the contribution from 
the cross-correlation term with blazars. In Sec. |IV[ 
we present the predicted angular power spectrum of 
DGRB from Galactic subhalos, and find combined con- 
straints on the annihilation cross section using extra- 
galactic and Galactic contributions. We conclude in Sec- 
tion |V] Throughout the paper, we adopt a flat cold dark 
matter model with a cosmological constant (ACDM) with 
the following cosmological parameters: fi™ = 0.277, 
flA = 0.723, Ho = lOO/ikms^^Mpc^i with h = 0.7, 
Us = 0.96, and ag = 0.81. 



II. EXTRAGALACTIC CONTRIBUTION 

In this section, we discuss the mean intensity and 
anisotropy of DGRB from dark matter annihilation in 
extragalactic halos. Much of the calculations are based 
on our earlier work [3 HH], but with extensions of the 
framework and significant updates on input models as 
explained below. 



A. Mean intensity 

The mean intensity of gamma rays from dark matter 
annihilation is given b}{^ 



liE)^ / dxW{[l + z]E,x){S^), 



(1) 



where E is the energy of photons, x is the comoving 
distance to a source at redshift z [x and z are used inter- 
changeably through the relation, dx/dz = c/H{z)], {S"^) 
is the variance of the overdensity field, S = {p— {p))/{p), 
and is often also referred to as an "intensity multiplier." 
W{E,z) is the window function that contains particle- 
physics information such as a velocity-averaged annihi- 
lation cross section times relative velocity, (crw), a dark 
matter mass, mdm, and a gamma-ray spectrum per an- 
nihilation, dNy^g^n^/dE: 
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where fidm = 0.23 is the density parameter of dark mat- 
ter, and pc is the critical density of the present Uni- 
verse. Here, t{E, z) is the optical depth for a gamma 
ray emitted at energy i5, for which we adopt the model 



^ Here we define the mean intensity, I, as the number of photons 
received per unit area, unit time, unit energy range, and unit 



solid angle, i.e., I{E) = dN/{dAdtdEdn). 



of Ref. [20] • Note that the annihilation cross section re- 
quired to produce dark matter at the right relic density 
by the thermal freeze-out mechanism with S-wave anni- 
hilation is {(Tv) ~ 3 X 10~^^ cm^ s""'^, which is largely 
independent of the dark matter mass [I]. 
The intensity multiplier is 



1 



^mPc 



dM 



X / dVpl,,,{T\M), 



(3) 



where M is the virial mass, dn/dM is the halo mass 
function, for which we adopt an ellipsoidal collapse 
model [31] [35], phost{r\M) is the density profile of a host 
halo of mass M, r is the comoving radius from the halo 
center, and bsh{M) is a boost factor due to annihilation 
in subhalos. 

We adopt an Navarro-Frenk- White (NFW) [33^ profile 
for host halos, 



Phostir) 



{r/rs)ir/rs 



1 



(4) 



where Ps and are the scale density and the scale radius, 
respectively, and this relation holds out to the virial ra- 
dius, Tvir, which, in turn, is given as a function of M and 
z through the relation: M = 4:nr^^j.Ayi,-{z)pc{z) /3, with 
Avir(z) = 187r2-h82(i- 39(^2 and d = r2,„(l-f z)3/[17„(l + 

zf + riA]-! [Mj. 

The scale radius is defined as Ts = T'vir/cvir, where 
Cvir(Af, z) is the concentration parameter, for which we 
adopt the model of Ref. [35] for masses below 2.5 x 
1O^'*M0 and that of Ref. [33] otherwise. By taking the 
volume integral of the density profile, phost {f) , out to r^ir 
and equating it to M, we obtain the scale density as 



Ps 



M 
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(5) 



The volume integral of the density squared has an ana- 
lytic form: 



1 - 
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(6) 



The gamma-ray intensity is further boosted by anni- 
hilation in subhalos, which is represented by the boost 
factor, bsh{M), for which we adopt a fitting formula 
based on results of recent numerical simulations [29]: 
6sh ~ llO(M2oo/lO^^M0)°-^^, where M200 is an enclosed 
mass within a radius r2oo in which the average density is 
200 times the critical density; there is a simple relation 
between M200 and the virial mass M [37]. This boost 
is realized if the subhalo mass function extends down to 
Earth-mass scales, Mgh.min = 10^^ Af©, which is a typical 
cutoff scale for the neutralino dark matter [38l - l4l] . We 
note, however, that the boost factor strongly depends on 
the minimum subhalo mass chosen, 6sh oc M, 



-0.2 



sh,min 



Given that a wide range of minimum subhalo mass is still 
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FIG. 1. Contribution to the intensity multiplier, {(5^) [Eq. ^] 
from different mass ranges at various redshifts. 



FIG. 2. Contributions to the intensity multiplier, (5^) 
[Eq. (|3|], from host halos (dashed) and host halos and sub- 
halos (solid) at 2 = 0. The solid line is the same as that in 
Fig.0 



allowed, as small as ^10 ^"^Mq [32], the annihilation rate 
may be boosted even further. 

Figure [T] shows the integrand of Eq. ([S]) as a function 
of the virial mass, M, for various redshifts, z. Namely, 
this shows what fraction of ((5^) is contributed from 
which mass ranges. As can be seen from this figure, 
the dominant contribution comes from cluster-size ha- 
los {M ^ IQ^'^ Mq) at z ^ and galaxy-size halos 
{M ^ 10^^ Mq) at z ~ 2. This is mainly because of 
the boost due to annihilation in subhalos. In Fig. [2] we 
compare the contribution from the host halos (dashed 
line) and that from the host halos and subhalos (solid 
line) at z = 0. If there are no gamma ray from subha- 
los, then the mean intensity would be dominated by the 
smallest dark matter halos. 

Figure [3] shows the intensity multiplier ((5^) as a func- 
tion of redshifts for the case with and without sub- 
halo contributions. Here we multiply {S"^) by (1-1- 
z)^ / y/ flm{l + z)^ + 57a to show contributions to the 
mean intensity from different redshift ranges. One can 
see that the presence of subhalos boosts the intensity by 
a factor of 10 at low redshifts, and by a factor of ^ 2 even 
at z = 5. The dotted lines are further multiplied by the 
absorption factor, e^'^, for observed energies oi E ^ 10, 
20, 50, and 100 GeV. There is little absorption for pho- 
tons received below 10 GeV, but this effect is significant 
for energies above tens of GeV and should be taken into 
account. 

Figure |4] shows the predicted mean intensity of DGRB 
from dark matter annihilation with mdm = 100 GeV, 
(av) = 3 X 10~^^ cm^^ s~^, and the bb annihilation chan- 



nel. This model gives the dark matter contribution that 
is as large as 10% of the mean intensity measured by 
Fermi-LAT gS] at ~ 10 GeV. This contribution is 
quite significant given that even the most dominant con- 
tributors known to date, i.e., unresolved blazars, con- 
tribute to DGRB at around the same level EH. 



B. Angular power spectrum 

The angular power spectrum at a given multipole, £, 
is given by 

CiiE) = J ^W^{[l + z]E,z)Ps2(^k^^,zy (7) 

Following conventions of recent publications (e.g., 
Ref. [37]), our definition of the angular power spectrum 
[Eq. ^] has the units of (cm'^ s'^ sr"! GeV''^)^ srQ 
which is referred to as the intensity angular power spec- 
trum. In order to obtain the fluctuation angular power 



^ When we compare theoretical predictions with the data, we must 
integrate a gamma-ray intensity in a given direction over energy 
within a given energy bin. We do this by replacing the window 
function, W{[1 + z]E,z), in Eqs. ([T} and ^ with the window 
function integrated over a given energy range. This gives Ci in 
units of (cm~'^ sr"-*^)^ sr. 
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FIG. 3. Intensity multiplier as a function of redshifts. The 
solid line shows the total intensity multiplier from host ha- 
los and subhalos, while the dashed line shows the intensity 
multiplier only from host halos. The dotted lines show the 
total intensity multiplier multiplied by the absorption factors 
at various energies {E — 100, 50, 20, and 10 GeV from left to 
right). 



FIG. 4. Predicted mean intensity of DGRB from dark mat- 
ter annihilation via the bb channel. The dark matter mass 
is mdm ~ 100 GeV, and the cross section is at its canoni- 
cal value, {av) = 3 x 10~^^ cm"^ s~^. The solid line shows 
the extragalactic contribution, while the dashed line shows 
the Galactic subhalo contribution (no smooth Galactic con- 
tribution is included). The Fermi-LAT data [13] as well as 
the contribution from unresolved blazars are also shown for 
comparison. 



spectrum that has units of sr and is adopted in the ear- 
Uer papers (e.g., Ref. [3]), one simply divides the intensity 
power spectrum by the mean-intensity squared, l'^{E). 

Here, Ps2{k,z) is the power spectrum of the overden- 
sity squared, S^, which can be divided into one- and two- 
halo terms [3]: 



2h/ 



(8) 



The one-halo term correlates two points in one identical 
halo, whereas the two-halo term does that in two distinct 
halos. Correspondingly, the angular power spectrum is 
also divided into two terms: 



ah/ 



(9) 



These two terms of the power spectrum can be explicitly 



written as 

P^Hk,z) = 



^711 Pc 



dM'-^\uik\M)\^ 



(l + 5sh(M)) / dVpl,,{r\M) 



^mPc 
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dM'^^^^^^u{k\M)b,{M,z) 



X (l + 6sh(M)) / dVpl^^MM) 



Plinik^z), 
(11) 

where Piin(fc, z) is the linear power spectrum of the mat- 
ter density field S, and bi (M, z) is the linear halo bias [32] . 
The power spectrum of S^, Ps^, depends on profiles of 
density squared in a halo of mass M, u(r|A/), where 
u{r\M) is normalized such that its volume integration be- 
comes unity; u(k\M) is the Fourier transform of u{r\M). 

Fourier transform of the density-squared profile, 
■u{k\M), is the sum of the density-squared profiles of the 
host halo and subhalos weighed by a fractional luminosity 
of each component: 



u{k\M) 



{thost(fc|M) + bsh(M)i2,h(fc|M) 
l + fosh(M) 



(12) 
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FIG. 5. Fourier transform of the density-squared profile, 
u{k\M), for various host-halo masses. The solid lines show 
the total (host halo and subhalos) profile, whereas the dotted 
lines show the host-halo profiles. Note that the total u{k\M) 
is normalized to unity at — > 0, and thus the dotted lines do 
not approach unity at fc — >■ but approaches f/(f -|- 6sii) [see 
Eq. (pi]. 



Here, we ignore a contribution from the cross term, 
2/3host('')Psh('"), which is important only when phost(?') ^ 
Pshir) at the same radius, r. Given that spatial distri- 
butions of the host halo and subhalo contributions are 
quite different (the host halo being important inside the 
scale radius and the subhalos being important outside), 
this approximation is very good. (See Appendix for the 
contribution of the cross term.) 

Fourier transform of the host halo profile, Uhost{'') c>c 
/^host('')' analytic form [S], but here we use an even 

simpler fitting formula. 



Uhost{k\M) 



1 



[l + a(fcr,)2/6]fc' 



(13) 



with a = 0.13 and b = 0.7, which is largely indepen- 
dent of Cvir ,20) . On the other hand, we obtain the 
density-squared profile of subhalos, UshirlM), by depro- 
jecting the surface brightness profiles of numerical simu- 
lations [MIES], assuming spherical symmetry as follows: 



-1 -3/2 



for r < r2oo, 



(14) 

with T] = 2.78. Note that the distribution of subhalos is 
typically more extended than the density profile of the 



FIG. 6. Relative contributions to the one-halo power spec- 
trum Pg2 at 2 = as a function of masses and wave numbers. 
The lines are normalized at M = 1 Mq. 



host halo. Its Fourier transform is 



sh(A:|M) = A 



dx- 



sin(Ka;) 
(a;2 -f 1/16)3/2 KX 
64 K cos K + sin k 



173/2 k(k2+772) 



(15) 



where k = fcr2oo and A « 0.64 is the normalization con- 
stant such that Msh(O) = 1. 

Figure [5] shows u{k\M) for various host-halo masses, 
M. When u{k) is close to unity (for small fc), the halo 
can be regarded as a point source. On the other hand, 
when u{k) deviates significantly from unity at a given 
wave number fc, the source extension cannot be ignored 
at that wave number, and the power spectrum (especially 
the one-halo term) is suppressed. Figure 5] shows that 
u{k) is larger for smaller host halos, which are less ex- 
tended. It also shows that the contributions from subha- 
los are more important for larger host halos. As the dis- 
tribution of subhalos is more extended than the density 
profile of the host halo, the subhalo contribution domi- 
nates at small k and the host-halo dominates at large k. 
This makes a hump at scales corresponding to the scale 
radius, r^. In other words, annihilation from the smooth 
host-halo component dominates inside the scale radius, 
where subhalos are tidally disrupted. 

Figure [6] shows the integrand of the one-halo power 
spectrum of (52, dP^i'ik, < M)/d\nM [Eq. at ^ = 

for various wave numbers, k. (Note that the lines are 
normalized at M = 1 AIq, and thus it shows relative 
contributions rather than absolute.) The bulk of the con- 
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FIG. 7. The same as Fig. [g] for fc = 1 Mpc ^, but in compar- 
ison with the host-halo contribution. 




Redshift z 



FIG. 8. One-halo contributions to the angular power spec- 
trum, Cl'', at £ = 10, 100, and 1000 as a function of redshifts. 



tributions come from large-mass halos, and halos smaller 
than a typical dwarf size (M < 10^ Mq) do not make any 
sizable contributions to any relevant ranges of k. This is 
particularly true for scales larger than a typical cluster 
size (i.e., k < 1 Mpc~^). For smaller scales, on the other 
hand, a relative importance of large-mass halos is smaller, 
as large-mass halos are more extended and thus the power 
from them is suppressed (as also shown in Fig. [s]). 

Figure [t] is the same as Fig. |6]for A: = 1 Mpc^^, but 
the host-halo contribution is also shown. We find that 
the impact of subhalos on the power spectrum is much 
greater than that on the mean intensity (as shown in 
Fig. [2]), and it can boost the power spectrum by almost 
four orders of magnitude at this particular wave number. 

After projected on the sky, the three-dimensional wave 
number, k, and the angular multipole, £, are related by 
k = £/x{z) for a given redshift, z [see Eq. Q]. Note, 
however, that this simple relation is valid only for small 
angular scales, £ ^ 1 (e.g., [9l[l0]), on which we mainly 
focus in this paper. 

Figures [S] and [9] show the contributions to the angular 
power spectra at £ = 10, 100, and 1000 from one-halo 
and two-halo terms, respectively, as a function of red- 
shifts. To calculate a contribution to Ci from a given 
z, we multiply Ps2 by some combination of functions of 
redshift [see the integrand of Eq. ([t]) and also the redshift 
dependence in Eq. (2])]. We find that lower multipoles are 
dominated by nearby sources: one-halo terms a,t £ = 10, 
100, and 1000 are dominated by sources at z ^ 0.002, 
0.02, and 0.2, respectively, whereas the dominant contri- 
butions to the two-halo term come from somewhat higher 



redshifts. 

One should not, however, include contributions from 
arbitrarily small redshifts in the integral of Eq. ([t]), as 
cosmic variance in such small redshifts is so large that 
taking the ensemble average (as we do here) no longer 
makes sense. In addition, when a source is sufficiently 
close, it should give enough gamma-ray fluxes to be iden- 
tified as an individual source which we can remove from 
the map. In the following discussion, we use three dif- 
ferent minimum redshifts in the integration of Eq. ([7|: 
Zmin = 0.001, 0.003, and 0.01. The r.m.s. overdensity 
within radii corresponding to these three redshifts are 
1.4, 0.74, and 0.28, respectively. Since none of these are 
much larger than one, we can argue that our results by 
setting these lower cutoffs are not subject to strong cos- 
mic variance. We also note that, as we show below, this 
choice does not make any significant difference for the 
multipoles we consider {£ > 150) in this paper. 

In Fig. [TOj we show both the one-halo and two-halo 
terms of the angular power spectrum, Cf, integrated 
over 5-10 GeV energy range, for three different val- 
ues of 0min = 0.001, 0.003, and 0.01. The particle 
physics parameters are (av) = 3 x 10~^^ cm^ s""'^ and 
"T-dm — 100 GeV, and we assume the bb annihilation 
channel. Note that scales as (crw)^. Taking smaller 
Zmin increases the power at large angular scales, in par- 
ticular for the one-halo term, because of the contribu- 
tions from closer, more extended halos. For the rest of 
this paper, we shall use z^in = 0.003 for definiteness. 

The thick solid line in Fig. [TO] shows the upper limits 
on the angular power spectrum [27\ (with the blazar con- 
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FIG. 9. The same as Fig.[8]but for the two-halo contributions 
to the angular power spectrum, Cf^. 



tribution subtracted between £ = 155 and 504 for 5- 
10.4 GeV, Ci55<£<504 < 8 X IQ-^o (cm-^ g-i sr-i)^ sr. 
The current upper limits in 5-10.4 GeV, whose energy 
region is most sensitive to dark matter particles with 
TOdm = 100 GeV annihilating into 6& (see Fig. [4|, are 
three orders of magnitude larger than the prediction 
with the canonical particle physics parameters. Recalling 
Ci oc (crw)^, we find an upper limit on the cross section of 
(f ^ 8 X 10^^^ cm^ s^^ for the 100-GeV dark matter 
annihilating into hh. 

The blazar-subtracted upper limits on the angular 
power spectrum are available in several energy bands, 
1.04-1.99 GeV, 1.99-5 GeV, 5-10.4 GeV, and 10.4- 
50 GeV [2H]. In Fig. [TT| we show the combined upper 



limits on the annihilation cross section for the hh channel, 
{av)i,i, as a function of the dark matter masses, Wdm, us- 
ing all the available data on the power spectrum. (Note 
that the model is still based only on the extragalactic 
contribution.) We calculate the upper limit on the cross 
section such that the predicted Ci dX £ = 155 (the low- 
est multipole at which the measurement is reported) is 
equal to the 2(T upper limit reported by Ref. [55]. The 
limits from each energy range are shown separately as 
the dotted lines. The combined limits shown as the solid 
line are simply the best of the four limits at a given dark 
matter mass. A more optimal analysis would improve 
these limits. 
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FIG. 10. Predicted angular power spectrum of DGRB from 
dark matter annihilation in extragalactic halos, integrated 
over the energy range between 5 GeV and 10 GeV. (This 
energy range is most sensitive to dark matter particles with 
Wdm = 100 GeV annihilating into hh; see Fig. H) The one- 
halo (solid) and two-halo (dotted) terms are shown separately, 
for three different minimum redshifts ^min = 0.003, 0.001, and 
0.01 (from top to bottom lines). The particle physics param- 
eters are {av) = 3 x era? and mdm = 100 GeV, and 
the hh annihilation channel is assumed. For comparison, the 
thick solid line shows the upper limit on the angular power 
spectrum in 155 < I. < 504 from the 22-month data of Fermi- 
LAT |27l with the blazar contribution subtracted 1281. 



III. CROSS CORRELATION BETWEEN DARK 
MATTER AND BLAZARS 



As astrophysical gamma-ray sources, such as blazars 
associated with supermassive black holes at the center of 
galaxies, reside in dark matter halos, there is a spatial 
correlation between gamma rays from dark matter an- 
nihilation and blazars. The angular power spectrum of 
dark matter annihilation including the cross correlation, 
which is equal to the total power spectrum minus the 
blazar power spectrum, is given by 
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Cp — c. 
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(16) 
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FIG. 11. Upper limits on the annihilation cross section into 
the bb final state as a function of the dark matter masses, 
mdm. The limits are obtained from the measurement of the 
angular power spectrum a.t £ = 155 [27] with the blazar con- 
tribution subtracted [28], and the theoretical model of dark 
matter annihilation from extragalactic halos. The dotted lines 
show the Hmits from different energy ranges: 1.04-1.99 GeV, 
1.99-5 GeV, 5-10.4 GeV, and 10.4-50 GeV from bottom to 
top at low-mass region. The solid line shows the combined 
limits, and the horizontal dashed line shows the canonical 
cross section for the thermal-freezeout scenario. 



where is the cross-power spectrum computed frorrj^ 

C^{E)=. j %Wb{[1 + z\E,z)Wa,^{[\ + z]E,z) 
J X 



X 



(17) 



where the subscripts "i?" and "dm" denote blazars and 
dark matter, respectively. Here, Wdm is the same as 
Eq. ([2|, while Wb is the window function for the mean 
intensity of blazars [similarly defined as Eq. ([T]) but re- 
placing ((5^) with 1]: 



Wb{[1 + z]E,z)=x' 



dC <pE{C,z)TEiC,z), 
(18) 

where C is the differential luminosity (i.e., luminosity per 
unit energy) at a given energy E, Te is the differential 



^ We include blazars in the cross-correlation, ignoring other 
sources of gamma rays. This may be justified, as blazars are 
so far known to be the most dominant extragalactic gamma-ray 
sources in the GeV energy regime. Other promising sources in- 
clude star- forming galaxies |21l I46H48| . 



FIG. 12. Predicted angular power spectra of DGRB in 5- 
10 GeV from dark matter annihilation only (dotted), dark 
matter-blazar cross correlation (dashed), and the sum of the 
two (solid). The particle-physics parameters are the same as 
those in Fig. |10[ For the blazar bias, we use = 1 (bottom 
dashed/solid) and 5 (top dashed/solid). 



flux, and <&£;(£, z) is the gamma-ray luminosity function 
of blazars. The upper limit of the integral, C{Fetz), 
is the luminosity giving the flux corresponding to the 
point-source sensitivity of Fermi-LAT at a given redshift 



z, for which we adopt 10 cm 



above 100 MeV. 



More formal definitions of these quantities can be found 
in Ref. [TO]. 

The three-dimensional cross power spectrum is 



P^{k,z) = {bB{z)) 



1 



^mPc 



dM^^{M,z)u{k\M) 



X 6i(M,z)(l-f6sh(M)) J dVpl,,,{r\M) 

X Plin{k,z) 

= {bB{z))^P^Hk,z)PUk,z), (19) 

where {bB{z)) is the blazar bias averaged over the lumi- 
nosity function and the flux as follows: 



(bBiz)) = 



/ dC<pE{C,z)J^E{C,z)bB{jC,z) 
JdC^EiC,z)TE{C,z) 



(20) 



with the luminosity-dependent bias, bB{C,z), and the 
same upper and lower limits of integration as those in 
Eq. 



Note that this power spectrum [Eq. ( 19 )] includes the 
two-halo term only. While the one-halo term, where dark 
matter annihilation happens in the same halo that hosts 



9 



a blazar, also exists, the previous study TU] shows that 
the one-halo term of the cross correlation is much smaller 
than the two-halo term; thus, we shall ignore the one-halo 
term of the cross-power spectrum. 

For the luminosity function, z), we adopt 

the luminosity-dependent density-evolution (LDDE) 
model [TT1|49] with the gamma-ray spectra assumed to be 
a power law with an index of 2.4, which is in agreement 
with the spectrum of resolved blazars as well as that of 
the DGRBJ^ Compared with the earlier study [TT] where 
the luminosity function was based on pre-Fermi data, we 
here adopt different values for parameters of the lumi- 
nosity function (k = lO""*, q = 3.5, and 71 = 1.05) such 
that the model reproduces the flux distribution of blazars 
resolved by Fermi [44] . 

In Fig. [4] we show the blazar contribution to the mean 
intensity. It is difficult to explain the DGRB intensity 
measured by Fermi with blazars alone, in agreement with 
the previous study [S] [52] . 

Figure [12] shows the angular power spectra from the 
cross correlation with (&b) = 1 and 5, for the energy 
range between 5 GeV and 10 GeV. The particle-physics 



parameters are the same as those in Fig. 10 We find that, 
if the annihilation cross section is around the canonical 
value required to produce dark matter particles at the 
right abundance by the thermal-freezeout mechanism, 
then the cross-correlation term cannot be ignored. This 
is particularly important when the bias of blazars is as 
high as 5. 



As the dark matter term (C^ 



Ih 



2h 



') and the 



cross-correlation term {C^ oc (crv)) scale with (av) dif- 
ferently, one may ask, "At which value of the annihilation 
cross section does the cross-correlation term become im- 
portant?" Figure [T3] shows Ce=i55 divided by the upper 
limit, C}^5<-^<-5Q4, for 5-10 GeV. The dashed Hne is with- 
out the cross correlation, while the solid lines are with the 
cross correlation with bias of 1 and 5. We find that the 
dark matter term dominates at {av) > 3x 10^^^ cm'^ s^^, 
whereas the cross-correlation term dominates at lower 
cross sections. 



IV. GALACTIC CONTRIBUTION 

Annihilation signals from subhalos in the Galactic halo 
containing the Milky Way are typically comparable to, 
or even greater than, the extragalactic contribution [T] [H] 
[Tni[ISllin]; thus, we must also take the Galactic subhalo 
contribution into account. We shall follow an analytic 
treatment presented in Ref. 20J. 

On the other hand, we do not include the contribution 
from a smooth density profile of the host halo of the 
Milky Way in our calculation. In the anisotropy analysis 



More elaborated spectra in combination with the luminosity 
function and the DGRB intensity are studied in Refs ,50. ,51 . 
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FIG. 13. Predicted angular power spectrum at ^ = 155 in 5- 
10 GeV divided by the current upper limit, C}^5<^<504, as a 
function of the annihilation cross section, {pv) . The solid lines 
include the dark-matter-blazar cross-correlation with (bs) = 
1 and 5, while the dashed line does not. The vertical dotted 
line shows the canonical cross section for thermal WIMPs. 



of the Fermi-LAT data [27], the low Galactic latitude 
region of |&| < 30 deg is masked. Using a density profile of 
the smooth Galactic component of Ref. [20] , we find that 
this mask brings the smooth contribution to the mean 
intensity down to about 10% of the subhalo contribution, 
as shown in Fig. [14] As the smooth component does not 
have much power on small angular scales, it can be safely 
ignored for ^ > 100 (but it can be comparable to the 
subhalo and extragalactic contributions on large angular 
scales, £< 10). 

The angle-averaged mean intensity from dark matter 
annihilation in Galactic subhalos can be computed from 



Ish = / dL 



(21) 

{L) dL 

where L is the gamma-ray luminosity of a subhalo, s is 
the line-of-sight coordinates, and dn^/dL is the angle- 
averaged luminosity function of subhalos. The lower limit 
of the line-of-sight integral, (L), corresponds to the flux 
sensitivity of Fermi-LAT, i.e., L = 47rs^Fsons- 

The subhalo luminosity, L, is related to the subhalo 
mass, M , via 



L = B. 



sh- 



{(TV)N. 



7,ann 



9m 2 



dV,^ ft\(rsh|M), (22) 



where rvir,MW is the virial radius of the Galactic halo, 
the subscripts "sh" denote subhalos, and B^h is a boost 
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FIG. 14. The same as Fig. |4] but for the Galactic smooth 
component in the high Galactic latitude region outside of |6| > 
30 deg (solid). The dashed line is the same as that in Fig. [4] 



FIG. 15. Predicted angular power spectra of DGRB in 5- 
10 GeV from dark matter annihilation in extragalactic halos 
(dotted), Galactic subhalos (dashed), and the sum of the two 
(sohd). 



factor due to the presence of substructure in subhalos 
(sub-subhalos) . 

We assume that the density profile of subhalos, psh, is 
well described by the NFW function [Eq. Q]. Then, the 
volume integral of the density-squared has the analytic 
form given by Eq. ([6]) with the concentration parameter, 
Cvir, replaced with Ccut that corresponds to the cutoff 
radius of subhalos, i.e., Ccut = ^ cut /fs- With this mass- 



luminosity relation [Eq. (22 1] and the subhalo mass func- 
tion dnsh/dM, one can compute the luminosity function. 
Most model inputs such as the subhalo mass function, 
spatial distribution, and mass-concentration relation are 
adopted from recent numerical simulations of the Galac- 
tic halo, Aquarius [S^. More details on how to apply 
these models to gamma-ray computations are described 
in Ref. [20]. 

The intensity angular power spectrum is 



sh 



1 



167r2 



dL 



d^ 



dL 



s 



(23) 

where Ush(fcj M) is the Fourier transform of the density- 
squared profile of the subhalo distribution, which is given 
by Eq. ( 13 ) if the density distribution of subhalos follows 
an NFW profile. Note that Eq. (p3| only includes "one- 
subhalo" term, where one corre 
identical subhalo. There is 



ates two points m one 
however, the two-subhalo 
term that correlates two points in two distinct subhalos, 
but this term is much smaller than the one-subhalo term 
at small angular scales J20]. 



from Galactic subhalos and extragalactic halos (but not 
including the cross correlation). We have used the canon- 
ical model of the Galactic subhalos given in Ref. |53j . 
which has the mass resolution of about 4 x 10^ Mq- We 
have extrapolated their result down to the Earth-mass 
scale (model Al of Ref. [20 ). The intensity power spec- 
trum is about the same for both the extragalactic and 
Galactic components, with the latter slightly larger in 
the angular scales constrained by Fermi-LAT. 



In Fig. 16 we show the limits on {av) from the Fermi- 



Figure 15 shows the predicted angular power spectra 



LAT data, taking into account both the extragalactic and 
Galactic terms. As expected, the limits from either alone 
are similar, and the combined limits improve by a factor 
of two. In particular, for low-mass dark matter particles, 
the combined limits are only a factor of three larger than 
the canonical cross section. The limits are weaker for 
larger masses. 

While our limits are not yet as stringent as those ob- 
tained from analyses of dwarf galaxies [SH [55| or galaxy 
clusters [IS] [5^ , where the canonical cross section is al- 
ready excluded for low-mass ('-^lO GeV) dark matter par- 
ticles, they are not so far away (i.e., only a factor of 
three to four worse). Also, our limits are derived in a 
completely different way: they are based on the diffuse 
emission rather than on individual objects, and they are 
based on anisotropy rather than on the mean intensity. 
It is certainly encouraging that the first limits using the 
DGRB anisotropy are already not so far away from the 
best limits. 
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FIG. 16. The same as Fig. |11| but for the hmits obtained 
from the Galactic subhalos (dashed); extragalactic halos (dot- 
dashed); and the sum of the two (solid). The dot-dashed line 
is the same as the solid line in Fig. |11| The dotted lines show 
the Galactic-subhalo limits from each of four energy bins. 
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FIG. 17. Density-squared profiles of a host halo (dotted); 
subhalos (solid) ; and a host-subhalo cross-term (dashed) . The 
mass of the host halo is A4 = 10^'' Mq and the redshift is 
z — 0. The virial radius is 1.2 Mpc and the scale radius is 
210 kpc. 



V. CONCLUSIONS 



In this paper, we have used the angular power spec- 
trum of DGRB recently detected in the 22-month data of 
Fermi-LAT [57] to place limits on the annihilation cross 
section of dark matter particles as a function of dark 
matter masses. As dark matter annihilation occurs in 
all cosmological halos and subhalos, our model includes 
all the contributing terms in the extragalactic halos, the 
Galactic subhalos, and the cross correlation between dark 
matter annihilation and blazars. The smooth Galactic 
component is predicted to be sub-dominant in the high 
Galactic region (|6| > 30 deg) and is ignored. 

We have revised our earlier model of the extragalac- 
tic contribution by including the results from recent 
numerical simulations of the subhalo distribution 
Combined with the model of the Galactic subhalos of 
Ref. [20j, we find that the Galactic and extragalac- 
tic contributions are comparable to each other. The 
cross correlation with blazars is important for annihi- 
lation cross sections smaller than the canonical value 
{{av) < 3 X 10^26 gj^3 g-i)^ 

By comparing our model with the upper limit on the 
non-blazar contribution to the angular power spectrum of 
DGRB [55], we find upper limits on the annihilation cross 
section as a function of dark matter masses as shown 
in Fig. |16| The current limit from anisotropy excludes 
regions of (av) > 10~^^ cm^ at the dark matter mass 
of 10 GeV, which is only a factor of three larger than the 



canonical value. The limits are weaker for larger dark 
matter masses. The first limits from DGRB anisotropy 
that we find in this paper are already competitive with 
the best limits in the literature. 

Our limits will improve as Fermi collects more data. 
At the same time, an improvement in the analysis can 
significantly improve our limits. Currently, the angular 
power spectrum on large angular scales, £ < 155, is not 
used because of a potential contamination by the Galac- 
tic foreground emission (such as pion decay). As the 
angular power spectrum of DGRB from dark matter an- 
nihilation, (without multiplying by rises towards 
low multipoles, including the low-multipole data will sig- 
nificantly improve the limits. This line of investigation 
(i.e., a better characterization and removal of the Galac- 
tic foreground) should be pursued. 
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FIG. 18. Fourier transform of the lines in Fig. |17[ 



Appendix A: Effect of host-subhalo cross-term 



The dark matter annihilation signal is proportional 
to the density squared. When both host halo and sub- 
halo contributions are present, one has the host-density- 
squared term, p^osti ^^e subhalo-density-squared term, 
Pgjj, and the host-subhalo cross term, 2phostPsh- In our 



analysis, we have ignored the cross term [see Eq. ( 12 )], as 



spatial distributions of the host halo and subhalo contri- 
butions are quite different (the host halo being important 
inside the scale radius and the subhalos being important 
outside). In this Appendix, we quantify the importance 
of the cross term. 



Figure [17] shows the density-squared profiles of a host 
halo, subhalos, and the cross term. As expected, the 
cross term becomes comparable to the other terms only 
within a narrow window in radii. For the host halo mass 
of M = 10^^ Mq and the redshift of z = 0, the cross term 
becomes comparable to the other terms at r '-^ 20 kpc. 



which is 1/10 of the scale radius, 
shows the Fourier transform. 



210 kpc. Fig. 18 
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